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We study the pion mass dependence of the p(770) and /o(600) masses and widths from one and 
two-loop unitarized Chiral Perturbation Theory. We show the consistency of one-loop calculations 
with lattice results for the M p , f n and the isospin 2 scattering length 020- Then, we develop and apply 
the modified Inverse Amplitude Method formalism for two-loop ChPT. In contrast to the /o(600), 
the p(770) is rather sensitive to the two-loop ChPT parameters - our main source of systematic 
uncertainty. We thus provide two- loop unitarized fits constrained by lattice information on M p , /„-, 
by the qq leading \/N c behavior of the p and by existing estimates of low energy constants. These 
fits yield relatively stable predictions up to m, ~ 300 — 350 MeV for the p coupling and width as 
well as for all the /o(600) parameters. We confirm, to two- loops, the weak jtv dependence of the 
p coupling and the KSRF relation, and the existence of two virtual /o(600) poles for sufficiently 
high m^. At two loops one of these poles becomes a bound state when m n is somewhat larger than 
300 MeV. 

PACS numbers: 



I. INTRODUCTION. 

The spectrum of the lightest meson resonances in 
QCD, particularly the scalars, is still not well understood 
from first principles. Lattice QCD can provide, in prin- 
ciple, a rigorous way to extract non-perturbative quanti- 
ties from QCD, but current calculations are hindered by 
difficulties, like those associated to disconnected graphs, 
and their need to use relatively high quark masses, so 
that appropriate extrapolation formulas are needed. In 
Ref. [l[ the quark mass dependence of the mass and width 
of the p and a (or / (600)) mesons was predicted. Both 
the mass and width were obtained from the poles in 7T7r 
scattering amplitudes generated with the one-loop In- 
verse Amplitude Method (IAM) [2M|. The IAM is ob- 
tained from a dispersion relation based on analyticity, 
elastic unitarity and Chiral Perturbation Theory (ChPT) 
Q. The relation between and the quark mass is 
well known and model independent within ChPT, which 
also provides the correct m w dependence of low energy 
scattering amplitudes up to a given chiral order. There 
were, however, two aspects in which the calculation in [I[ 
needed further improvement: on the one hand, at that 
time we were only able to compare with a single lattice 
p(770) mass calculation at one value of m^. On the other 
hand, we presented a calculation up to one-loop, namely 
next to leading order (NLO) unitarized ChPT, so that 
the main source of systematic uncertainty, coming from 
higher orders, could not be estimated. 

In this paper we address these two issues. First, in the 
next subsection we briefly review the notation and the 
one-loop unitarization formalism used in 1]. Next, in 
Sect [Til we show the good agreement of the one- loop pre- 
vious calculation with several lattice group results, not 
only for the p(770) mass at many values of the pion mass, 
but also for and the isospin-2 scattering length pion 
mass dependence. In addition, since ref.|f| provided a 
strong indication that the pirn coupling was almost quark 



mass independent, we perform here an explicit calcula- 
tion of the p and a couplings from the residues of their 
associated poles, confirming the very weak dependence of 
the p pole on the pion mass, in sharp contrast with the a 
behavior. Next, in Sect lllll we present our calculations 
for unitarized two-loop Chiral Perturbation Theory. Un- 
fortunately, the dispersive formalism used in IfJ, which 
is a modification of the well known IAM 0, Q that in- 
corporates correctly the Adler zeros, was developed in 
[J] in detail only for the one-loop case. Thus, in subsec- 
tion HIO] we provide the justification and the expression 
that modifies the IAM to incorporate the Adler zeros 
correctly up to two loops. Next, in subsection IIII Bl we 
discuss the large uncertainties in the two-loop low en- 
ergy constants, that will dominate our systematic errors, 
particularly for the p(770), as explained in subsection 
IIII CI Therefore, to ensure the correct mass dependence, 
in subsection IIII Dl we fit these constants not only to ex- 
perimental data but also to the existing lattice results 
with the additional constraint of respecting the qq lead- 
ing order 1/N C behavior for the p as well as existing esti- 
mates of the low energy constants. We show that all these 
constraints can be accommodated with fairly reasonable 
low energy parameters. Finally, in subsection IIII El we 
provide predictions for the pion mass dependence of the 
controversial /o(600) scalar resonance parameters, which 
are remarkably stable under the two-loop uncertainties. 
In addition, we provide predictions for the coupling con- 
stant of the p and the Kawarabayashi-Suzuki-Riazuddin- 
Fayyazuddin (KSRF) relation. 



A. IAM and mlAM 

The p and a resonances appear as poles in 7T7t scat- 
tering partial waves with definite isospin / and angular 
momentum </, in the / = 1, J — 1 and I = 0, J = 
channels respectively. Elastic unitarity implies for these 
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partial waves and physical values of s: 

Im£(s) = cr(s)\t(s)\ 2 => Im l/t(s) = -a(s), (1) 

where s is the Mandelstam variable and a(s) = 2p/y/s, p 
being the center of mass momentum. Consequently, the 
imaginary part of the inverse amplitude is known exactly. 
ChPT amplitudes, being an expansion £ = £2 + £4 + • • • , 
with tk — 0(p k ), satisfy unitarity only perturbatively: 

lm£ 2 = 0, Im£ 4 = cr|£ 2 | 2 , ■■■■ (2) 

Let us recall that £2 corresponds to a tree level calculation 
with the leading order (LO) chiral Lagrangian, whereas 
£4 contains the one loop diagrams with LO vertices, plus 
tree level terms from the NLO chiral Lagrangian. The 
LO Lagrangian has no free parameters, but just m v and 
/„•. Higher order Lagrangians contain Low Energy Con- 
stants (LECs) that renormalize the loop divergences and 
whose values contain the information about the underly- 
ing theory, QCD. These LECs carry a scale dependence 
to cancel that of the loop integrals, so that observables 
are scale independent and finite order by order. 

The IAM [J] uses elastic unitarity and the ChPT ex- 
pansion to evaluate a once subtracted dispersion rela- 
tion for the inverse amplitude. The analytic structure of 
l/£ consists on a right cut (RC) from threshold to 00, 
a left cut (LC) from —00 to 0, and possible poles com- 
ing from zeros of £. The scalar waves vanish at the so 
called Adler zero, 5,4, that lies on the real axis below 
threshold. Its position can be approximated with ChPT, 
•5a = s% + S4 + • • • , where £2 vanishes at s 2 , £2 + £4 van- 
ishes at s 2 + S4, and so on. We can write then a once sub- 
tracted dispersion relation for 1 /£, the subtraction point 
being s A , 

1 = s~s A r Iml/£(s') 
£(s) 7T J RC 3 («' - s A )(s> - s) (3) 
+ LC(l/t) + PC(l/t), 

where LC(l/t) stands for a similar integral over the left 
cut and PC(l/£) is the contribution of the pole at the 
Adler zero. Note that, as l/£ already has a pole at s A 
the usual subtraction constant terms are actually part of 
the pole contribution term. 

On the right cut we can evaluate exactly Im 1/t = 
—a = — Im£4/£ 2 , as can be read from eqs. ([I]) and 
([2]). Since the left cut is weighted at low energies we 
can use ChPT to approximate LC(l/t) ~ LC(— t±/t\). 
The pole contribution PC (1/t) can be safely calculated 
with ChPT since it involves derivatives of £ evaluated at 
s A , which is a low energy point where ChPT is perfectly 
justified. Altogether, we arrive to a modified one-loop 
IAM (mlAM) formula 0]: 

■mlAM _ L _2 

£ 2 - £4 + A mIAM ' 

A mIAM = w S; \ _ p2 ~ S A )(S - S 2 )[t' 2 (s 2 ) - £ 4 Q 2 )] 

s - S A 

(4) 



where the prime denotes a derivative with respect to s 
and where we use for s A in the numerical calculations its 
NLO approximation s 2 + S4. The standard IAM formula 
is recovered for A mIAM — 0, which is indeed the case 
for all partial waves except the scalar ones. In the origi- 
nal IAM derivation @, Q A mIAM was neglected since it 
formally yields a higher order contribution and is numer- 
ically very small except near the Adler zero. However, 
if A mIAM is neglected, the IAM Adler zero occurs at s 2 , 
correctly only to LO, is a double zero instead of a sim- 
ple one, and a spurious pole appears close to the Adler 
zero. All of these caveats disappear with the mlAM, and 
the differences between the IAM and the mlAM in the 
physical and resonance region are less than 1%. 

It is important to remark that, in the above deriva- 
tion, ChPT has not been used at all for calculations of 
t(s) for positive energies above threshold. Note that the 
use of ChPT is well justified to calculate s A and PC, 
since these are low energy calculations. ChPT has also 
been used to calculate the left cut integral, which, despite 
extending to infinity, is heavily weighted at low energies, 
which once again justifies the use of ChPT. The approx- 
imation of the left cut and the subtraction constants up 
to a given order in ChPT - with the subtraction point 
chosen at low energy -, together with the elastic approx- 
imation are the only approximations used to derive the 
IAM from the dispersion relation, but no other model de- 
pendent assumptions are made. In particular there are 
no spurious parameters included in the IAM derivation, 
but just the ChPT LECs, and f„. 

Remarkably, the simple Eqs.@ (either the IAM or the 
mlAM) ensure elastic unitarity, match ChPT amplitudes 
at low energies, and, using LECs compatible with existing 
determinations, describe fairly well data up to somewhat 
less than 1 GeV, generating the p, K* , a and k resonances 
as poles on the second Riemann sheet. 

Of course, other unitarization techniques are possible, 
but in order to improve the NNLO IAM these would im- 
ply the use of coupled channels, which are not needed 
for the a and the p, higher orders of ChPT, which, if 
available, could be incorporated in a higher order IAM 
version, or a different approximation of the left cut. For 
the latter, we can distinguish two possible regions for im- 
provements: either at low energies, where we could sys- 
tematically include more orders of ChPT, again leading 
to a higher order IAM version, or at high energies (which 
are nevertheless suppressed by the subtractions) where 
the left cut should be modeled introducing more - non 
ChPT - parameters whose dependence on QCD is not 
known, or further assumptions beyond those used in the 
IAM. There is also, of course, an ambiguity on where to 
choose the subtraction point, but, as shown in differ- 
ent choices contribute to higher order corrections in the 
subtraction constants and the pole contribution terms, 
which numerically differ very little from one another. 

Finally, there are other unitarization techniques which 
are very successful and simpler than the IAM. However 
these can be recast, in the elastic regime, as the IAM 
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plus further a ppr oximations, like dropping crossed and 
tadpole terms [26[ - and therefore keeping a spurious pa- 
rameter like a cutoff or an unknown subtraction constant 
to regulate the theory or keeping just the leading order 
Q , in which case the p cannot be reproduced with natu- 
ral size parameters (see discussion below). Some of these 
simpler methods can be easily understood in terms of 
resummations of particular sets of Feynmann diagrams, 
but without the said additional simplifications there is no 
proof that the IAM, and much less so the mlAM, can be 
obtained from a simple diagrammatic resummation. All 
these simpler methods are known to provide physical re- 
sults for the scalars rather similar to the IAM. However, 
in order to extrapolate to non-physical quark masses one 
would need additional model dependent assumptions on 
the behavior of the spurious parameters. For the rea- 
sons explained in the last two paragraphs the IAM is 
the most adequate and complete technique to study the 
quark mass dependence of the p and a elastic resonances. 

For our purposes in sec lIII Al it is important to remark 
that it has been shown [5[ that the scalars can actually 
be generated mimicking the LEC, tadpole and crossed 
channel diagrams by a cutoff of natural size, and thus it 
is said that scalars are "dynamically generated" from, es- 
sentially, meson-meson dynamics (meson loops). In con- 
trast, to generate the vectors, a precise knowledge of the 
LECs is needed, namely, of the underlying, non meson- 
meson QCD dynamics. As we will see this makes the 
p(770) much more sensitive to the still poorly known two- 
loop LECs, whereas the sigma is rather stable under such 
higher order corrections. 

Thus, by changing m w in the IAM amplitudes we can 
study how the generated p and a poles evolve, so we 
can predict the dependence of their masses, widths and 
couplings on m w . The values of to be considered 
should lie in the applicability region of ChPT and should 
allow for some tttt elastic regime. Both criteria would 
fail above, roughly, ~ 500 MeV, which was taken as 
the upper bound for one- loop calculations in [l[ . We will 
see here that the approach is not reliable much before, 
namely, around m T ~ 300 — 350 MeV at least. 

In Q the mlAM was used for the p and a chiral ex- 
trapolation, because, for the scalar and at high m^, one 
resonance pole gets near the IAM spurious pole, a prob- 
lem that is nicely solved with the ml AM. Nevertheless, 
in the physical region and near the other generated poles, 
the differences between IAM and mlAM approaches are 
almost negligible, even for high pion masses. What we 
have briefly reviewed is just the NLO, or one-loop, case 
derived in Q. For this work we will need the two-loop 
version of the mlAM that we will explicitly calculate in 
sect IIIII below. But first we will show that the predictions 
obtained in [l| for the chiral extrapolation of the one-loop 
case are in quite good agreement with recent lattice re- 
sults on the pion mass dependence of the p mass, /„• and 
isospin 2 scattering length. 



II. ONE-LOOP RESULTS 

Within the SU(2) ChPT formalism, for tttt scattering 
only four LECs appear at 0(p 4 ), and are denoted by 
Since for now we just want to compare the 
predictions in [l[ with lattice results, we will use the 
same values as in [lj. Namely, we use lO 3 ^ = 0.8 ± 3.8, 
lOHl = 6.2 ± 5.7 from [|; and l\ and l r 2 are obtained 
from a mlAM fit to phase shift data up to the resonance 
region, 10 3 Z[ = -3.7 ± 0.2, 10 3 ^ = 5.0 ± 0.4. All the 
LECs are evaluated at fi = 770 MeV. 

Let us first compare with lattice results and then we 
will calculate explicitly the coupling constant of reso- 
nances to two pions. 



A. Comparison with lattice 

Around the time of the publication in [l| , several lattice 
calculations were published providing results for the pion 
decay constant f„ @, [H, [lj] and also for the tttt scalar 
isospin 2 scattering length 020 [13] • What we will show 
next is that, despite these lattice results were not used 
as input in the calculations of they are fairly well 
described within our one-loop formalism. 

In particular, when changing in our amplitudes we 
have to change accordingly the value of f„ , as it also de- 
pends on m n . Note that, of course, the f„ calculation 
is not unitarized, but just standard ChPT. However, / ff 
is an important factor in all our unitarized calculations. 
Actually, since we are using ChPT up to 0(p 4 ) in the 
IAM dispersion relations, we evaluate the f v dependence 
on to that order. At 0(p 4 ), the only LEC that ap- 
pears in the dependence of / T is l\, whose value is 
fixed here to that given in as commented above. In 
the top panel of Fig. [1] we compare the resulting one 
loop dependence of f n on with some lattice results. 
The grey area covers the uncertainty in l\ only, and one 
should also recall that lattice uncertainty bars are sta- 
tistical. We consider the spread of the different lattice 
groups as an estimate of their systematic errors. With 
these remarks in mind, we can see that the depen- 
dence implemented in our approach is compatible with 
that calculated from the lattice. Our grey band does 
not extend beyond = 0.5 GeV because, at the very 
least, that is for sure an applicability bound for our cal- 
culations. In the next sections, we will see that at those 
pion masses the two-loop uncertainties are actually too 
big to make any significant quantitative claim with our 
method, which we will find reliable only up to roughly 
m w ~ 300 - 350 MeV, at most. 

The I — 2, J = tttt partial wave can be calculated 
with the IAM, and indeed it was included in the fit to 
phase shift data 1], thus constraining its energy depen- 
dence. However, the m n dependence of 0,20 was not con- 
strained with any input. Of course, it can be easily pre- 
dicted and is interesting to compare it with the available 
lattice data as a consistency check of the method. Thus, 
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FIG. 1: Top: One loop ChPT and lattice results for the f„ 
dependence on m^, but with the l\ parameter used in our 
IAM fits. Lattice points are from % 111 [3. Bottom: IAM 
and lattice results for the I = 2, J = scattering length 020 
dependence on m^. Lattice points are from [l2l 



in the bottom panel of Fig. Q] we show the lattice results 
for 020 compared to our IAM calculation. Once again, 
one should take into account that the mlAM error band 
only covers the uncertainties in the one-loop LECs. In 
view of the figure, and taking into account that our curve 
is not a fit to these data, we consider that our predictions 
are in fairly good agreement with these lattice results. 

Of course, we could refit our approach including the 
data in FigfTJ and we would be getting a better agree- 
ment, but at this point we only want to check the con- 
sistency of our results and that we do not need a fine 
tuning to describe lattice results, which will be included 
as input of two-loop fits in the next sections below. 

We have thus checked our method for consistency 
against available lattice results on two quantities other 
than resonance masses, which we address now. First, 
let us recall that the mass M and width T of a narrow 
resonance are related to its pole position as y/s po i e = 
M — iT/2, and this notation is usually also kept for wide 



FIG. 2: Comparison of our IAM results for the M p (calculated 
here as the point where the phase crosses tt/2) dependence 
on m n with some recent lattice results 0-Qj]]. The grey band 
covers only the error coming from the LECs uncertainties. 



resonances, as done in However, lattice calculations 
do not provide results in the complex plane. For that rea- 
son, in this work we will also consider for the p(770) the 
most usual and physically intuitive definition: namely, 
that the mass of the resonance corresponds to the energy 
where the scattering phase shift reaches tt/2. This is 
the value where the modulus of the scattering amplitude 
shows a peak, and we will thus call it "peak mass" . Of 
course, this definition is only valid for narrow resonances, 
and is a very good approximation for the p, which is the 
one for which more reliable lattice results exist. For the 
a we will stick to the pole mass definition. 

Thus, in Fig.[5]we show the results of our IAM calcula- 
tion of the p mass evolution when varies, which is dis- 
played as a grey band that covers the uncertainties in the 
one-loop fitted LECs only. This behavior agrees nicely 
with the estimations for the two first coefficients of the 
M p chiral expansion [15| . This figure is relatively similar 
to Fig. 3 in |l| , except that we plot versus and not 
and that we have defined the mass as the energy where 
the phase shift crosses tt/2. This latter choice ensures 
that the physical p mass from the PDG [39| lies within 
our uncertainty band. In contrast, the "pole" mass would 
lie somewhat lower. This difference between "pole" and 
"peak" masses decreases as m w increases, and they are 
almost indistinguishable around 350 MeV. Of course, we 
are now comparing with a compilation of lattice results 
from different collaborations. Let us remark that due to 
the finite lattice volume, for some of the lattice results, 
the minimum energy with which pions are produced is 
larger than the resulting M p and therefore the resonance 
has zero width. With these caveats in mind, and in view 
of the large systematic deviations between different lat- 
tice collaborations, one can conclude that our one-loop 
prediction shows a rather good agreement with the bulk 
of lattice results. 
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FIG. 3: p (continuous line) and a (dashed line) couplings 
to 7T7T evolution with m^. Since the two conjugate a poles 
become two real poles after reaching the real axis, there is 
a coupling for the pole that moves towards threshold ("up") 
and another one for the pole moving away ("down"). 

B. Resonance couplings 

In Ref. [l| it was shown that both the p and er widths 
- calculated with the mlAM from the imaginary parts of 
the pole positions- decrease as is increased. It was 
also shown that the decrease in T p is largely kinematical, 
following remarkably well the expected reduction from 
phase space as M p approaches threshold, whereas the a 
width decreases in a quite different way of that provided 
only by phase space reduction. Following this argument 
it was concluded in [l| that the effective coupling of the 
p to 7T7T must be almost independent and that the a 
coupling should show a strong dependence. However, 
we did not perform an explicit calculation of these cou- 
plings, and we will provide it here. For elastic amplitudes 
with a given isospin / and angular momentum J, we can 
define their coupling g to two pions as 

9 2 ~ lim ( S ~ S f ^ S \ (5 ) 

s^s p -p ZJ 

where s p is the position of the resonance pole on the 
second Riemann sheet, and p is the center of mass mo- 
mentum. The above definition, if used as such, is very 
unstable numerically. For that reason we have calculated 
the residue applying Cauchy's Theorem to a small circle 
surrounding the pole. 

In Fig. [3] we show the results for the couplings. The 
continuous line shows the p7T7r effective coupling evolu- 
tion as m-K changes, defined as in Eq.© from the IAM 
partial wave and normalized to its physical value. As ex- 
pected, it is almost independent, since it deviates by 
less than 5% from its original value when changing the 
pion mass from 139.57 MeV up to 450 MeV. 

The results of our explicit calculation thus confirm 
quantitatively the p coupling independence suggested 



qualitatively in 1]. This result is very relevant because it 
justifies the constancy assumption made in lattice studies 
of the p(770) width 

In contrast, the <77T7r coupling, plotted as the dotted 
line of Fig. [3] shows a strong dependence. The dra- 
matic peak in this curve, around 320 MeV corresponds 
to the pion mass where, as shown in the two conju- 
gated poles of the sigma meson on the second Riemann 
sheet join on the real axis below threshold. As the pion 
mass increases beyond that value, the two poles remain 
real and we therefore have to plot two different branches. 
One of these two poles (labeled "lower" ) stays away from 
threshold, at least within the mass range of this study. In 
contrast, the other one (labeled "upper") moves towards 
threshold and, eventually jumps onto the first Riemann 
sheet. For the one loop calculation this occurs slightly 
above m w = 450 MeV, but we will see that for the two- 
loop calculations this may occur for pion masses even 
only slightly above 300 MeV. Let us also remark that it 
is a nice consistency check to see that, as the pole reaches 
threshold, the coupling tends to zero, as it is expected 
from general arguments [H, . 

As it could be expected, our results both for the /?(770) 
and the a are also in very good agreement with a very 
recent one-loop calculation using unitarized SU(3) ChPT 
\{\ , since the strange quark mass plays a very little role 
in the p(770) and a values of their masses and widths. Of 
course, we expect our SU(2) results to be more reliable 
than those of SU(3) since, as it is well known, the SU(2) 
ChPT convergence is much better. The small differences 
with respect to |l6j can only be attributed to small dif- 
ferences in the choice of LECs. Actually the SU(2) LECs 
that we use here were fitted in [l| to pion-pion scattering 
data only. In contrast, in [l6[ the SU(3) LECS used there 
were fitted to the whole elastic meson-meson scattering 
data below 1 GeV plus some lattice results. The effect 
of the the kaon loops, which have been integrated out in 
the large kaon mass limit, is negligible compared to our 
uncertainties. 



III. TWO-LOOP FORMALISM AND RESULTS 

In this section we extend our analysis to the next to 
next to leading order (NNLO), namely, two- loop unita- 
rized ChPT. As we will see, the present knowledge about 
the two-loop LECs is rather poor. However, it will be 
enough to obtain sufficiently robust predictions for some 
observables, particularly those related to the a meson, 
and to estimate the size of the uncertainties due to higher 
order terms in the different parts of the calculation. First 
of all, of course, we have to rederive the mlAM formalism 
for the two-loop case. 
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A. The two-loop mlAM formalism 

As we have already commented, instead of the IAM, 
we need the mlAM in order to describe properly the sub- 
threshold region in the scalar waves. As explained in 
subsection II Al in the dispersion integral of Eq. (J3j) we 
now have to evaluate the left cut and the pole contribu- 
tion using the 0(p 6 ) expansion of the ChPT amplitudes: 
t = t 2 + <4 + te + .... Of course, now the position of the 
Adler zero has to be evaluated also up to that order, i.e., 
sa — s 2 ± S4 + sq. After a tedious calculation along the 
lines of (H, we arrive at the 0(p 6 ) version of the ml AM: 



A' 



rlAM 



tIAM 



t 2 -U + t 2 /t 2 -t 6 + A mIAM ' 

2t4(s 2 )t 4 (s2) tl{s 2 ) 



U{s 2 ) 



+ t 6 {s 2 ) + 



t' 2 (S2) 

(s - s 2 )(s A - s 2 



SA 



t' 2 (s 2 )(s - s 2 ) 
- ( t' 2 (s 2 ) - t' 4 (s 2 ) 



t'±{s 2 ) 2 + t'l{s 2 )U{s 2 ) 
t' 2 (s 2 ) 



(6) 

As a technical remark, let us note that we have to in- 
troduce an additional subtraction to ensure convergence. 
Note also that we will now need the second derivative of 
£4 with respect to s at s 2 . Again, the standard two-loop 
IAM 0, [HI is recovered for A mIAM = 0, which occurs 
for all waves with J > 0. We can now use the 0(p 6 ) 
mlAM to calculate the chiral extrapolation of the p and 
a resonance poles, and check the consistency with the 
0(p 4 ) results. 



B. Two-loop ChPT and the low energy constants 

The two-loop 7T7T SU(2) ChPT scattering amplitude 
was calculated in [24[ and contains six additional LECs. 
These are denoted r\ , ■ • • r§ and their values are poorly 
known. Moreover, this 0(p 6 ) calculation is expanded in 
terms of (m^ / fn) 2 , where f n is the pion decay constant 
evaluated at the physical value of the pion mass, not in 
the chiral limit. That is fair enough to describe the phys- 
ical scattering amplitude, but in order to extract the pion 
mass dependence of the scattering amplitude one also has 
to include the 0(p 6 ) ChPT /„. pion mass dependence [29j : 



Jo Jo 



k 

fo 



1 

"N 



L [711 + 2l 2 ~ 11 



29 



12NJ 



it 
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-211H 



1 

N 2 



13 
192 



(7) 



where L = -k log (to 2 //i 2 ), N = I671" 2 and f/ is the rel- 
evant combination of 0(p 6 ) LECs that appears in f^, 



TABLE I: Sample of LECs: First row: Roy-Steiner eqs. SU(3) 
analysis of nK scattering. Second and third rows: K14 anal- 
ysis to 0(p 4 ) and 0(p 6 ), respectively. Naively, we have com- 
bined quadratically the SU(3) LECs errors there. Fourth row: 
Roy Eqs. analysis with uncertainties from imaginary parts 
and unknown 0(p 6 ) LECs combined quadratically. Fifth row: 
Roy Eq. analysis of low-energy irir scattering up to two loops, 
whose errors "only account for the noise seen in their calcula- 
tions". Last row, values used in [l[ with the IAM. All LECs 
are evaluated at the scale fi — 770 MeV 



Analysis 


10 3 ll 


w 3 r 2 


io 3 zs 


io 3 z2 


0(P 4 ) [22J 
0(p 4 ) [30] 
0(p 6 ) [30J 
0(p 6 ) [20] 
0( P 6 ) f21] 


-4.9 ±0.6 

-4.5 
-3.3 ±2.5 
-4.0 ±2.1 
-4.0 ±0.6 


5.2 ±0.1 
5.9 

2.8 ± 1.1 
1.6 ± 1.0 

1.9 ±0.2 


2.1 
1.2 ± 1.7 

0.8 ±3.8 


17 ± 10 

5.7 
3.5 ±0.6 

6.2 ± 1.3 


IAM [1] 


-3.7 ±0.2 


5.0 ±0.4 


0.8 ±3.8 


6.2 ±5.7 



poorly known once again. For that reason, we will use 
lattice data on f v to stabilize it in our fits. 

Let us now remark that the two-loop leading log con- 
tributions, which are numerically dominant at low en- 
ergies, do not depend on the ri constants, but just on 
the one-loop l\. For this reason it is well known that 
the values of the l\ ■ ■ ■ l\ 0(p 4 ) LECS can vary sizably 
between the 0(p 4 ) and 0(p 6 ) analysis, as it is shown in 
Table [R where we quote the LECs obtained in several 
works |20l - [22l [30j . Let us remark that error bars in the 
l\- ■ -l\ usually correspond to the statistical uncertainty 
("noise") in the input, but there are other large system- 
atic sources of uncertainty that are most likely dominant. 
Hence, it should not be surprising that some of the 0(p 6 ) 
LECs deviate somewhat from the estimates when taking 
into account leading log terms from 0(p s ) and higher 
orders, as we will do with the IAM. 

Finally, it has been shown [25( that the LECs can be 
understood as the effective couplings that result from in- 
tegrating out heavy fields. Indeed, most one-loop LECs 
values are saturated by vector resonance exchange, like 
the p(770), but note that the a plays a little role, if any, 
in the actual values of the LECs. These resonance sat- 
uration estimates have also been extended to the 0(p 6 ) 
LECs [24j, but these are very uncertain and are custom- 
arily assigned a 100% uncertainty. We have collected 
them in Table UH together with some other estimates 
coming from dispersive analysis of pion-pion scattering 
& 

At this point we want to make clear that the LECs do 
not depend on the quark mass. This is a rather obvious 
statement for people familiar with ChPT but, when pre- 
senting in conferences and workshops [3ll . [32j previous 
results from Ref. [l[ or partial preliminary results from 
this work, we have found that people get confused since 
we have just stated that the l\ are saturated by reso- 
nance exchange - like p exchange - but we have shown 
that the p mass depends on the quark mass. The reason 
for this confusion is that resonance saturation is usually 
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interpreted as a ~ \jM R contribution to the l\, with Mr 
the physical mass of a resonance. However it is actu- 
ally a ~ X/M R0 contribution, with Mro the resonance 
mass in the chiral limit, namely M R = M R0 + 0(m%). 
Numerically, using Mr instead of Mrq makes a small dif- 
ference - neglected when obtaining LECs estimations - 
but is incorrect in terms of ChPT. A term like ~ V-^fl 
coming from integrating out a heavy resonance should 
be re-expanded as ~ 1/M| (1 + O(ml/M R0 )). The first 
term contributes to Z[, but the next 0(rn^) term counts 
as a higher order in ChPT and therefore does not con- 
tribute to the same order as k does. The same occurs to 
all orders, so, as stated above, the LECs do not have to 
be readjusted when the mass of higher resonances change 
with m v . They do not depend on the quark mass. 



C. Resonance sensitivity to LECs 

After reviewing briefly the standard ChPT two-loop 
calculation, we can now use it to generate the p(770) and 
/o(600) resonances, for which we use the mlAM. Thus, 
we show once again in Table Q] the LECs we obtained in 
the one- loop IAM fit in [l], which are fairly compatible 
with the non-unitarized determinations, lying roughly in 
between the one and the two-loop bulk determinations. 
Naively, this could correspond to the fact that the one- 
loop IAM, reproduces not only the one-loop ChPT ex- 
pansion but also the numerically relevant s-channel two- 
loop diagrams. 

At this point it is important to recall the different role 
that LECs play in the generation of the p(770) and the 
a resonances. In particular, it has been shown that the 
a can be easily generated within the chiral unitary ap- 
proach [f| from the leading order ChPT - which only de- 
pends on f n and m w - and a natural cutoff, whereas that 
is not feasible for th e p( 770), which needs the input from 
the one-loop LECs [2g]. This is also understood from 
the 1 /N c behavior of these resonances [27], |28[ : namely, 
the p(770) behaves nicely as a qq state when generated 
from the IAM - a N c dependence due to the the 1/N C 
behavior of the leading LECs. In contrast, the a does 
not behave predominantly as a qq and this dependence is 
mostly due to logarithmic terms, which are independent 
of the LECs and are generated from meson loops. This 
explains why the a, despite being lighter than the p, is 
not contributing so sizably to the value of the low energy 
constants. In summary, the terms containing LECs play 
a crucial role in the generation and location of the p pole 
in the IAM, but not so much for the a. This is due to the 
fact that to generate the p we need the input from the 
underlying QCD dynamics of a qq state, and we expect 
this to occur to all orders in ChPT, whereas the sigma is 
dominated by the scale f v [j| and should depend much 
more mildly on the underlying QCD dynamics encoded 
in the LECs. 

This is actually what we find when we look for the m w 
dependence of the p pole with the 0(p 6 ) IAM; its be- 



havior is quite unstable under the large uncertainties of 
the 0(p 6 ) LECs. If we leave completely free the 0(p 6 ) 
parameters within their huge estimated uncertainties, we 
cannot make any two- loop prediction for the p(770). For 
this reason, at two loops we will fit not only experimen- 
tal data but also the LECs values. In addition, by fit- 
ting the experimental data only, as was done in [lj, one 
constrains mostly the combinations of 0(p 6 ) chiral pa- 
rameters that govern the energy dependence. However, 
one cannot expect to constrain the 0(p 6 ) LECs combi- 
nations that govern the pion mass dependence that is 
of interest here. Fortunately, as we have seen in previ- 
ous sections, there is a large amount of lattice data on 
the pion mass dependence of the p(770) mass, the 1 = 2 
scattering length and /„■, which can be used to constrain 
further the 0(p 6 ) LECs and will be included in our fits. 
Once this is done, we will obtain predictions for the 
dependence of the p coupling and width, as well as on all 
the /o(600) parameters, which is where the most inter- 
esting discussion is still going on. 



D. Unitarized two-loop constrained fits 

We have thus fitted the mlAM to the elastic pion-pion 
scattering phase shifts shown on the left column of FigJU 
Let us remark that we have fitted data up to 1 GeV 
for the (1,1) and (2,0) waves and up to 800 MeV for 
the (0,0) channel; beyond that energy the effects of the 
/o(980) are important and cannot be reproduced with 
our single channel formalism. There are, of course, cou- 
pled channel unitarization formalisms [H, H3| , which are 
very successful and generate the / (980) among other 
resonances, but they lie beyond the scope of this work, 
mostly because of their simpler treatment of the left cut 
- if dealing with it at all-, or their dependence on ad- 
ditional parameters, all of which can introduce further 
model dependences. 

In addition, in order to constrain further the LECs that 
govern the pion mass dependence, we have also fitted the 
lattice results shown in Figs.l and 2 in previous sections, 
although only up to a pion mass of 350 MeV. 

Still, we have 11 parameters to fit, and even with these 
constraints, because of the large correlations between pa- 
rameters, during the fitting procedure some LECs can 
take values very far from their typical ones, even of dif- 
ferent order of magnitude, for tiny improvements in the 
data \ 2 . For that reason, we have also considered an 
averaged Xleg's t erm as a constraint for our fits, that 
measures how far the LECs are from some reference val- 
ues that we provide in Table [TTJ To further constrain 
the fit, we will also require the p mass and width leading 
1/N C scaling to follow a qq pattern, so we also consider 
a X 2 p ~qq measure, as described in [28|], to constrain the 
p behavior to that of a qq. Note that uncertainties for 
N c = 3 are of the order of 30%, but this constraint be- 
comes stricter as N c grows. Nevertheless, we will never 
apply this constraint for N c larger than 20, since oth- 
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erwise the theory would become weakly interacting and 
the whole unitarization procedure would loose sense, as 
repeatedly explained in [3ll HH . 

In summary, we are considering several % 2 -like func- 
tions that, when smaller or close to 1, ensure a good de- 
scription of each feature described above: x 2 a ta> XM°** fce > 
XIT 1 "' X 2 f l attlce X 2 p -- qq and xIecs- The problem is' that 
many of these are not really well defined x 2 functions in 
the statistical sense. The reasons are that some uncer- 
tainties we use are theoretical (as for the 1/N C behavior), 
and that some sets, both for real data or lattice, are in- 
compatible with each other and we have to guess some 
systematic uncertainty of the order of the difference be- 
tween different sets. In addition, the number of "data 
points" to be fitted for each feature is very different, and 
we could get a bad description of one feature with few 
"data" at the expense of a tiny improvement on another 
feature with more "data" , but affected by systematic er- 
rors crudely estimated. Thus, there is not a single fit 
of data and lattice minimizing the sum of all x 2 func- 
tions, since we do not know how to weight each one of 
them against the others. We are nevertheless presenting 
four different "fits" A, B, C and D where we have im- 
posed that each one of the x 2 should be relatively close 
or smaller than 1. This is still quite a strong constraint, 
and ensures, as we will see in Fig. 21 that all features are 
fairly well described up to the applicability region. Thus, 
in Table IPTl we show the parameters of four different fits, 
which, as seen in FigHJ cover the data and lattice results 
fitted up to pion masses of 350 MeV. In particular, we 
show the M p dependence on versus the lattice data, 
the I J = 00, 11, 20 partial waves, as well as the 0,20 and f„ 
dependence on m T compared with lattice results. Note 
that these fits tend to prefer the somewhat stronger Mp 
pion mass dependence found by the MILC 0], ETMC [|] 
and RBC-UKQCD [§] collaborations. The leading 1/N C 
behavior of the p pole is also shown in Fig. [5] 



Fit 


Fit A Fit B Fit C Fit D 


Reference 


io a zr 


-5.0 


-4.7 


-5.0 


-4.0 


-3.3 ±2.0 


10 3 Z£ 


1.7 


0.95 


1.7 


1.24 


1.9 ± 1.0 


iq% 


0.82 


0.82 


-6.0 


0.82 


0.82 ±3.8 




6.5 


4.96 


3.5 


6.5 


6.2 ±2.0 


10 4 n 


-0.6 


-1.0 


-0.7 


-0.6 


-0.6 


10 4 r 2 


1.3 


1.3 


3.7 


1.5 


1.3 


10 4 r 3 


-1.7 


-0.29 


2.7 


-3.3 


-1.7 


ioV 


2.0 


4.2 


2.8 


0.95 


-1.0 


10 4 r 5 


2.0 


2.3 


2.0 


1.7 


1.5 


10 4 r 6 


-0.56 


-0.98 


-0.5 


-0.7 


0.4 


ff 


-3.4 


-1.8 


-2.3 


-4.6 






TABLE II: LECs for the 0(p 6 ) fits A, B, C and D and the 
reference values we use in XlbCs • F° r the 0(p 4 ) LECs we have 
used values that cover the different sets in Table I, whereas 
for the 0(p 6 ) LECs we show estimates from [2l|, [24|] 

In view of the different x 2 P er data points (abbrevi- 
ated as x 2 /(#P ints)) that we list in Table ITU and will 
describe in detail below, we consider that fits A and D 



x7(#points) 


Fit A Fit B Fit C Fit D 


S 00 Sol. B [25J 


0.8 


1.1 


0.5 


0.6 


<5 00 Sol. C [35] 


5.6 


1.1 


2.6 


5.7 


5 00 below 400 MeV 


0.87 


0.8 


0.9 


1.0 


5u 


1.2 


1.6 


0.7 


0.7 


$20 


0.2 


1.0 


0.3 


0.4 


IVlp 


0.77 


0.68 


0.88 


1.35 


^lattice 
"20 


0.2 


4.0 


0.06 


0.2 


rlattice 
JTT 


1.4 


1.1 


0.6 


0.8 


N c 


1.4 


0.6 


1.2 


1.5 


LECs 


1.5 


3.8 


3.2 


1.4 



TABLE III: x 2 /(#P°i n ts) for each feature fitted. See main 
text for details of the 'x calculation. 



are the best ones, particularly because their LECs are 
quite compatible with the values of the 0(p 4 ) LECs and 
reasonably close to the crude expectations for the 0(p 6 ) 
ones. The difference between them is that fit A describes 
better the M p pion mass dependence - even up to very 
high pion masses - and not so well the 020 scattering 
length, whereas fit D does the opposite. In addition, 
we provide two other fits to illustrate the uncertainties. 
First, on fit B we try to force the softer M p pion mass 
dependence found by the CP-PACS- Adelaide [lOj] and 
QCFSF [ll| collaborations. Note that this fit then prefers 
the less negative phase shift data of the (I, J) = (2, 0) 
channel, an, as a consequence it does not describe very 
well the mass dependence of the 020 scattering length. 
Also, the agreement with the LECs estimates is worse 
than for fit A, as seen in Table [III In addition we arc 
showing the results for fit C, which is very similar to A 
and D up to pion masses of 300 MeV, but shows a rather 
dramatic drop of the p mass above 350 MeV. It is a gen- 
eral feature of our fits that they can be made reasonably 
compatible with data up to 350 MeV, but beyond that 
point they can start diverging widely for not very large 
changes in the LECs. That is the reason why we consider 
that our approach is not reliable beyond a pion mass of 
300 to 350 MeV, depending on the observable. 

Before proceeding to the next section we want to pro- 
vide, for each feature included in our fit, the x 2 P er num- 
ber of data points. Of course, in view of the clear incom- 
patibilities of different sets of data for a given quantity - 
from experiment or lattice -, we have either added some 
systematic uncertainty to cover different sets, or fitted 
some particular subsets of data, as we detail next. Let 
us emphasize that, since we use the I AM and have made 
several approximations on its derivation, we are obviously 
not aiming at precision. As it is seen in Figf4]the spread 
of our four fits covers roughly the different data points. 

The chaotic situation with the (0,0) phase shifts above 
450 MeV (obtained from irN — > irirN experiments) is 
well known and very visible in FigHJ Even the same 
CERN/Munich experiment [35| provides 5 different anal- 
ysis (A, B, C, D and E) incompatible with each other. 
Among them, data sets B and C have been shown to be 
the ones that satisfy better several dispersive constraints 
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FIG. 4: Different fits to physical nix scattering phase shifts Su and lattice results on the p mass, / T and the 020 scattering 
length. For the lattice results, we have not extended our fits to the grey area {m^ > 350 MeV, although it is displayed to show 
how the fits deteriorate, or not, beyond 350 MeV. The phase shift data comes from [33T - |35| ] and the lattice results from 0-0 ■ 



36]. Thus, we have added linearly 2 degrees of system- 
atic uncertainty to the points of these two sets - not much 
given the huge incompatibilities. In Table [TTTI we provide 
the resulting x 2 / (#points) for each fit with respect to the 
data sets B and C in [35[ • Note that, for each fit, at least 
one of the two experimental sets is well described, i.e., 
with x 2 /(#points) < 1. If we wanted x 2 /(#points) ~ 1 



for all data sets plotted in Fig|4]simultaneously, we would 
have to add 5 degrees as a systematic error, which again 
would not be much taking into account that differences 
between data sets in the 700 to 800 MeV region can be as 
high as 15 de gree s. Fortunately, the K14 decay data be- 
low 400 MeV [34( is of much better quality, but is easily 
reproduced by our fits without the need to add system- 
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FIG. 5: N c dependence of p(770) mass and width for the 
different fits described in the text. 

atic uncertainties as seen also in Table Mil 

The fits to other quantities are less complicated. De- 
spite they look relatively close, the two sets of (/, J) = 
(1, 1) 7T7T scattering phase shifts are not compatible, since 
they provide tiny statistical errors only. Their difference 
is of the order of two degrees, which once again we have 
added linearly as a systematic error. Concerning the 
(2, 0) wave, once more it is clear that there are incompat- 
ible sets of data, and we have added the same 2 degrees of 
systematic uncertainty. The resulting x 2 /(#points) for 
all our fits are given in Table Hill Note that, once again, 
fit B can be used to differentiate the (2, 0) "low data" set 
from the "high data" set and that it does not reproduce 
so well as the others the physical p shape. 

Concerning the lattice predictions for the p mass, there 
is also an obvious conflict between different collabora- 
tions, with differences in the M p value as high as 70 
MeV. Hence, the x 2 /(#points) given in Table Hill has 
been calculated for the points with below 350 MeV, 
adding a systematic uncertainty of 35 MeV. Note that 
if we took into account data up to = 400 MeV, the 
X 2 /(#P°ints) would have been 0.95, 0.93, 0.99 and 2.1 
for fits A to D, respectively. 

In the case of the /„■ lattice results, different collabora- 
tions have points clustered sufficiently close in groups of 
three to estimate a systematic uncertainty as half the dif- 
ference between the highest and lowest data point, which 
is added linearly to the statistical errors. The result- 
ing \ 2 for each fit with this prescription is shown in Ta- 
blelnH Note that we find some difficulty in describing the 
JLQCD results simultaneously with the other features. 

The leading order 1/N C behavior of the p has been 
adjusted to be exactly that of a pure qq state. This is 
believed to be a very good approximation, although given 
its large physical width it may easily have some small tttt 
component, that we neglect. Using the % 2 definitions and 
estimating the uncertainty to be exactly 1/N C times the 
leading term, as explained in (28[, we find the averaged 
X 2 for the p as a pure qq state shown in Table IIIII 



The last quantity that we have fitted is the (2,0) scat- 
tering length, where we find one data point right at 
niir = 350 MeV and four other points below. There are 
two collaborations and their points are relatively consis- 
tent with each other. Thus, and despite the large dif- 
ferences between different collaborations for other ob- 
servables, we have not added any systematic error to 
these points. For the four points strictly below 350 
MeV we show the x 2 /(#pomts) for each fit in table HTTl 
Once again fit B has serious problems describing this 
observable, but note that it describes better the "high 
data" set of (2,0) phase shifts. If we now include the 
point at tci-k = 350 MeV with its tiny uncertainty, the 
X 2 /(#points) of all fits except fit D grow beyond 4.5. 

It is clear that our fits A, C and D give a very good 
general description of, at least, one set of data for each 
observable up to 300 MeV, but not so good up to 350 
MeV, where they start deviating from each other also 
for the M p quark mass dependence. For that reason we 
will consider our IAM approach to be valid only up to, 
roughly, the 300 to 350 MeV region. This is, of course a 
crude estimate and varies from one observable to another. 

Finally, concerning the Xlecs> ^ us remark that to 
calculate it we use the reference values given in Table 
HI! Since the values of the 0(p 6 ) parameters come 
from resonance saturation estimates (with resonances of 
angular momentum smaller than 2), we have assumed 
that they are only correct in the order of magnitude and 
therefore have a 100% uncertainty. For all fits, most of 
the deviations come from the 0(p 6 ) LECs, which are the 
worst known. Note that, in general, but particularly in 
this case, fits A and D are the ones with better x 2 an d 
that is the reason why we consider them as our best fits. 
Fits B and C are just given for illustration of different 
scenarios and to show that the predictions that we will 
detail next are robust even allowing for larger uncertain- 
ties or deviations from our best fits. 



E. Predictions from the fits 

Once we have obtained a relatively good description 
of the data and existing lattice results on certain observ- 
ables, we will now use these three fits to obtain predic- 
tions for other observables. In general, the spread be- 
tween our curves should be considered as a naive indica- 
tion of our systematic uncertainties. 



1. The p width and coupling 

First of all, in Table IIVI we provide the values of the 
p pole position, s p , in the lower half plane of the second 
Riemann sheet and its coupling to two pions defined from 
the (1,1) partial wave as 

3 

9 P = -16?^ lirr^ ( s - s p) h\{s)j^, ( 8 ) 
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FIG. 6: Predictions for the pion mass dependence of p pa- 
rameters from the IAM two-loop fits described in the text 
and Fig|3] Top: We show the rho width from each IAM fit 
together with the behavior expected only from phase-space 
variation due to the changing p mass. Middle: the p coupling. 
Bottom: the Mpj % ratio. The last two are remarkably in- 
dependent of the quark or pion mass, which ensures that the 
KSRF relation is not spoiled within the applicability region 
of our approach. 





^fs~p (MeV) 


1 1 

\9p\ 


fit A 


754 - i 74 


6.1 


fit B 


772 - i 71 


5.9 


fit C 


759 - i 72 


6.0 


fit D 


763 - i 71 


5.9 



TABLE IV: Values of the p pole position in the lower half 
plane of the second Riemann sheet and pmr couplings for the 
different fits described in the text. 



where the normalization factors are chosen to recover the 
usual expression for the two-meson width of the p: 



r P = \g P f 



6vr A/2 



(9) 



If we approximate M p = Re J~Sp~ in the above equations, 



and we use the value of the coupling in Table IIV1 we find 
that the width of the p for fits A to D are 149, 144, 145 
and 142 MeV, respectively, to be compared with the PDG 
value of 149 ± 1 MeV. [39] . A rather good description, 
given the approximation, the quality of the data, and 
that we only provide central values. 

On the top panel of FigJS]we show the evolution of the p 
width for the two-loop unitarized ChPT fits described in 
the text and FigJU Although hard to see because they are 
almost overlapping, we provide, together with the result 
of each fit, the expected variation from phase space due 
to the change in M p only, assuming a constant p — 2ir 
coupling in a Breit-Wigner form. As it happened in the 
one-loop case [l|, we see that this constancy assumption 
of the p coupling is a very good approximation, so that 
this feature is rather robust under higher order ChPT 
corrections when fitting to data and lattice results on /„■ 
and a,2o- 

Moreover, in the middle panel of Fig|51 we show the 
actual two-loop IAM calculation of the coupling from the 
residues of the amplitude at the resonance pole. We can 
see that it is rather independent of the pion mass up to 
run ~ 350 MeV. Namely, fit B barely changes at all, fits 
A and D change by roughly 5% and only fit C changes 
by 15%, to be compared with the factor of 2.5 times that 
the pion mass is increased from its physical value up to 
350 MeV. Actually, this corresponds to an increase in the 
quark mass larger than a factor of 6. Thus, the relatively 
weak dependence of the p—2ir coupling on the quark mass 
is confirmed at two loops. This is of particular relevance 
to some lattice calculations that have assumed a constant 
coupling in studies of the p width [17]. Nevertheless, 
this prediction of M p with the quark mass is limited to 
pion masses below roughly 350 MeV, and more uncertain 
within our two-loop approach than at one-loop, as we see 
from the spread between different fits. Actually, for the 
less favored fit C we see up to a 15% variation, but note 
that the M p mass in this fit C starts behaving rather 
weird around 350 MeV, and most likely it is, very roughly 
speaking, only reliable up to the region between 300 to 
350 MeV, depending on the observable. 

Finally for the p(770), we show in the lower panel of 



12 





^s CT (MeVJ 


|.g CT (MeVJ 


fit A 






fit B 


474 _ 7 948 


O.tJ 


fit C 


Af\(\ _ 7 94^ 


o.o 


fit D 


AVX — i 971 


3.5 
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R. Garci'a-Martm et al. [40] 


474±6-i(254±4) 


3.58 ±0.03 


R. Kaminski et al. [42] 


442 - i290 


2.47 ±0.45 


J.A. Oiler [43] 


(443 ± 2) - i(216 ± 4) 


2.97 ± 0.04 



TABLE V: Values of the a pole position in the lower half 
plane of the second Riemann sheet and amr couplings for the 
different fits described in the text and some references in the 
literature. 



Fig|5]the evolution of the M p / f„ ratio, which comes al- 
most independent of m OT in all of our two-loop fits. This 
ratio is not really a prediction, since we have fitted both 
M p and /„• to lattice data. However, it is of particular 
interest for the well known KSRF relation [Hi, which 



provides an striking connection between the p - 
pling and the M p j ' f K ratio: 



9% 



M 2 p /8fl 



2tt cou- 



(10) 



and holds fairly well for the physical values of these con- 
stants. Since we have just checked the almost constancy 
of 

9p-rrm then, an almost constant M p / ratio as found 
in Fig|5] means that the KSRF relation also holds rather 
nicely, at least within our applicability region. This cor- 
roborates the one-loop results already found in unitarized 
SU(3) ChPT 0. 



2. The a parameters 

Let us turn to the predictions for the pion mass depen- 
dence of the a or /o(600) scalar resonance, whose nature 
is still very controversial. As we did with the p we pro- 
vide first the values of the resulting a pole positions and 
couplings in Table [V] obtained as follows: 
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-16-7T lim (s — s a ) too(s). 



(11) 



For comparison we also provide in the last rows of Ta- 
ble [V] results from other works in the literature. Let us 
remark that, without using the recent and very precise 
K14 decay data [34| . the a pole was obtained more than 
fourteen years ago at ^fs^ = 440 — i 245 MeV, with the 
single channel one loop IAM in |[. 

The 0(p e ) results are in quantitative agreement with 
the 0(p 4 ) ones for pion masses lower than about 300 
MeV. For instance, as the quark mass is increased, the 
relative growth of the a mass, defined as the real part of 
the a pole position, is slower than the pion mass growth, 
but still somewhat faster than for the p mass. 

As we saw for the one-loop case in previous sections 
and in [l|, as the pion mass grows, since the sigma mass 
grows slower, its width becomes narrower and narrower, 



and its two conjugate poles approach the real axis. But 
contrary to the p case, it is possible to have poles below 
threshold on the second Riemann sheet, and the a con- 
jugate poles actually meet at some point on the real axis 
below threshold. If the pion mass keeps on increasing, 
both poles stay on the real axis, but one moves very lit- 
tle, remaining at masses lower than threshold, whereas 
the other one increases its mass, until it eventually jumps 
to the first Riemann sheet. All these features occur once 
again at two loops and this double real pole structure is 
nicely seen in the top panel of FigJT] as a double branch 
for each fit. 

Let us remark that the appearance of two branches is 
not an artifact of the IAM, but is a general feature of 
scattering theory of scalar amplitudes with poles close to 
threshold [49| , also seen in other contexts |50j . It is just 
the way scalar poles approach to threshold as one changes 
the features of the interaction. Namely, there are no re- 
strictions on where a scalar pole should be on the real 
axis below threshold on the second Riemann sheet, ex- 
cept that poles appear in conjugate pairs out of the real 
axis, or on the real axis below threshold. In the first 
case, they obviously have the same "mass", but if this 
pair reaches the real axis, the two poles no longer have 
to be conjugated and hence the double branch is a gen- 
eral feature. In contrast, all non-scalar waves have cen- 
trifugal p 2J factors relevant around threshold, that force 
their second sheet poles to reach the real axis precisely 
at threshold [4!| where one of them jumps into the first 
sheet whereas the other stays in the second, as it happens 
here with the p(770). 

The IAM, of course, not only reproduces this general 
feature, but also provides an estimate of the pion mass 
where this apparent splitting occurs, which is not generic, 
but a specific value due to the QCD dynamics underlying 
the properties of the lightest scalar meson. 

Let us remark that all our fits are very consistent with 
each other for the er, despite their differences for the p 
behavior. As explained above, this is due to the fact that 
in the generation of the a the chiral loops play a domi- 
nant role. At NLO these are independent of the LECs, 
and at NNLO they only depend on the 0(p 4 ) LECs, but 
not on the . Therefore the a avoids most of the largest 
uncertainties that affect the p and, as a consequence, 
the results for the a are much more robust. There are 
however, some quantitative differences with the one-loop 
results in [l[: for instance, the point where the two con- 
jugate poles of the a meet on the real axis occurs for 
masses about 20 MeV lower; namely, at m, = 280 — 310 
MeV for the different two-loop fits versus 300 — 330 MeV 
for the one-loop description. This is a rather small correc- 
tion to the one-loop result and confirms the robustness of 
our results for the <r, even under higher loop corrections, 
at least up to ~ 300 — 350 MeV, depending on the 
observable. However, for higher pion masses the quanti- 
tative spread is much larger, although the four fits yield 
the same qualitative predictions for the poles in the two 
branches. Closely related to the decrease of the "split- 



13 




200 300 400 

m„ (MeV) 



500 




200 300 
m„ (MeV) 



500 



FIG. 7: Predictions for the a from the IAM two-loop fits 
described in the text and Fig(4] Note they are remarkably 
stable despite the differences between fits for the p. We show 
as a grey area the region where > 350 MeV, which marks 
the applicability bound of the approach. Top: we show the 
predicted pion mass dependence for the a pole mass. Note the 
appearance of two branches around 300 MeV, as explained in 
the text. Bottom: we show the dependence of the a coupling 
to two pions, which, as explained in the text, is rather strong 
particularly as the pion mass approaches the value where the 
two conjugate poles on the second Riemann sheet reach the 
real axis. For each fit, the thick line corresponds to the upper 
branch and the thin one to the lower branch. 



sion relations for the inverse amplitude, so that this pole 
on the first Riemann sheet is a zero for the inverse ampli- 
tude and therefore does not alter the analytic structure) . 
Our two-loop results seem to indicate that a conventional 
bound state - not a virtual one- might be found for pion 
masses higher than 290 — 350 MeV, contrary to the 460 
MeV we found at one- loop, which as we have seen was 
for sure outside the region where our approach is reliable. 
This is in qualitative agreement with some recent lattice 
results in [45], where they seem to find a bound state 
for m w ~ 325 MeV. Let us nevertheless recall the caveats 
raised from the very authors of (45| , since they cannot cal- 
culate accurately the width, and some possibly relevant 
contributions - mainly the disconnected contractions- 
have not been included in the calculation. Other lat- 
tice studies [46| have also suggested the existence of a 
"tetraquark" component for m, ~ 180 — 300 MeV. Let 
us note, however, that the binding energy of the states 
we find at two-loops seems to grow faster with the pion 
mass than for the one-loop case and on the lattice. How- 
ever, this occurs already in the region > 350 MeV, 
where we do not consider our approach reliable and the 
uncertainties are huge as seen by the spread of the fits in 
FigEl Unfortunately, other relevant lattice calculations 
for the a [13, 5H , lie beyond our reach. 

Finally, in the lower panel of FigJTJwe show our results 
for the a — 2tt coupling, obtained from the residue of 
the second Riemann sheet pole - or poles when there 
are two branches. The qualitative behavior is similar to 
the one-loop case shown in Fig[3l with a dramatic rise 
up to a peak that occurs at the pion mass where the 
two conjugate poles meet on the real axis. From that 
value onwards we thus have to draw two branches for 
each fit, and it can be noticed that the coupling for one 
of these branches reaches zero. This corresponds to the 
pion mass where the upper branch pole reaches the 27T 
threshold. The fact that at threshold the coupling goes 
to zero is in good agreement with the well known result in 
[HI]. Actually, this can be checked analytically because, 
as shown in |19j, the coupling is inversely proportional 
to the energy derivative of the one-loop function which 
is divergent at threshold. 



ting point" is the fact that the pole of the " upper branch" 
reaches the threshold faster than in the one-loop case as 
the pion mass grows. Note that the threshold variation 
corresponds to the line labeled 2m n in the figure, that 
the upper branches of all fits touch very soon after the 
two branch splitting, around ~ 290 — 350 MeV, versus 
TOtt ~ 460 MeV for the one loop calculation in [l(. 

The relevance of these results is that, when the up- 
per branch pole reaches threshold, it jumps into the first 
Riemann sheet, and becomes a usual bound state. (One 
might wonder if the dispersion relation for the IAM ap- 
plies now that there is a pole on the first Riemann sheet, 
but note that the IAM derivation is obtained from disper- 



IV. SUMMARY 

Using the IAM, which is based on analyticity, elastic 
unitarity and ChPT, we generate the poles associated 
to the p and a resonances without any assumption on 
their existence or nature. The IAM implements the pion 
mass dependence of observables through the subtraction 
constants up to a given order in ChPT. Thus, we can 
predict the dependence of the p and a pole positions on 
TOtt , as done in [l[ . Here we present new results that were 
missed in the previous paper. 

First, using the one-loop formalism, we have made a 
comparison of our previous results with some recent lat- 
tice data, showing that they are compatible. We have 
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also calculated the m, dependence of the pmr and aim ef- 
fective couplings, calculated from the pole residues, find- 
ing that the pirn coupling is almost m, independent, 
whereas the ct7T7t coupling shows a strong m n depen- 
dence. 

Finally, we have extended to two-loops the modified 
Inverse Amplitude Method formalism to account prop- 
erly for Adler zeros, which has been applied then to the 
0(p 6 ) calculation. Although no robust predictions can 
be made for the p mass, mostly due to the large uncer- 
tainties in the low energy constants, we have been able 
to describe the elastic scattering phase shift data and 
lattice results on f n and 020 with several fits with fairly 
reasonable values for such low energy constants, and the 
correct qq leading 1 /N c behavior of the p. 

With these fits we have obtained relatively robust pre- 
dictions for other p observables and all a parameters, at 
least up to ~ 300 — 350 MeV. In particular, we have 
confirmed the relatively weak dependence of the p — 2n 
coupling and the approximate validity of the KSRF re- 
lation. Concerning the sigma, whose results are much 
more robust than for the p since it has a much weaker 



dependence on the ChPT low energy constants, we have 
confirmed the appearance of two virtual poles for suffi- 
ciently high pion masses. One of these poles becomes a 
virtual state for m„ between, roughly, 300 and 350 MeV. 
We hope these results could be of use as a guideline for 
future extrapolations of lattice results down to physical 
quark mass values. 
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